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' We study theoretically the Josephson effect in junctions based on unconventional superconduc- 

^ , tors with diffusive barriers, using the quasiclassical Green's function formalism. Generalized bound- 

O ' ary conditions at junction interfaces applicable to unconventional superconductors are derived by 

' calculating a matrix current within the circuit transport theory. Applying these boundary condi- 

tions, we have calculated the Josephson current in structures with various pairing symmetries. A 
^ : number of predictions are made: (a) nonmonotonic temperature dependence in d-wave supercon- 

■ ductor/diffusive normal metal/d-wave superconductor (D/DN/D) junctions, (b) anomalous current- 
pH I phase relation in p-wave superconductor/diffusive normal metal/p-wave superconductor (P/DN/P) 
Q , junctions, (c) second harmonics in D/DN/D and P/DN/P junctions, (d) a double peak structure of 
^ ' the critical current in D/DF/D junctions, (e) enhanced Josephson current by the exchange field in 

I ' S/DF/P junctions. We have also investigated peculiarities of the Josephson coupling in D/DF/D, 

P/DF/P and S/DF/P junctions. An oscillatory behavior of the supercurrent and the second harmon- 

■ ics in the current- phase relation are studied as a function of the length of the diffusive ferromagnet. 

c/3 ; 

. 

■ I. INTRODUCTION 
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The Josephson effect^ has been studied in various types of junctions^i^i^ motivated by fundamental interest and 
(— I ^ potential applications for future technology. In superconductor / diffusive normal metal / superconductor (S/DN/S) 
Q junctions the critical current increases monotonically with decreasing temperature^i^i^i^ because proximity effect is en- 
O , hanced at low temperatures. Superconducting junctions with ferromagnetic interlayers have shown rich physics due to 
the interplay of proximity effect and the exchange field^'^. When DN is replaced by a diffusive ferromagnet (DF), it was 
' predicted that tt junctions can be realized^ i^°i^^i^^i^^i^"^i^^i^^ . The physical reason for a 7r-state is nonzero momentum 
J> , of induced Cooper pairs in the ferromagnelii^, similar to the so-called Fulde-Ferrel-Larkin-Ovchinnikov statei^ii^ in a 
magnetic superconductor. SFS tt junctions were realized experimentally by several group a-'^^'^^i^-'^'^^i^^'^'^i^^i^^'^'^i^^'^^ . 
CO , In d-wave superconductor junctions, one of the most remarkable phenomena is the formation of midgap Andreev res- 
' onant states (MARS) at interfaces^. The MARS stem from sign change of pair potentials of d-wave superconductors^. 
^ I ' In d-wave superconductor / insulator / d-wave superconductor junctions, 7r-junctions emerge due to the formation of 
' the MARSSS,!^. In order to clarify the role of proximity effect and MARS, Tanaka et al. have extended the circuit 
theory'^^ to the junctions with unconventional superconductor a^^i'^^'^^ . The conservation of matrix current enables one 
"Y^i ■ to apply the generalized Kirchhoff's rules to unconventional superconducting junctions and to derive the boundary 
' conditions for the Usadel equation^* widely used in diffusive superconducting junctions. Application of this theory to 
^ , the DN/d-wave superconductor (DN/D) junctions has shown that the formation of MARS strongly competes with the 
I • proximity effect in D N^^i^^ . It was also demonstrated that the formation of MARS coexists with the proximity effect 
' O _ in DN/p-wave superconductor (DN/P) junctions, which produces a giant zero bias conductance peak (ZBCP)'^^. 
^ ' Recently this theory has been extended to diffusive Josephson junctions with unconventional superconductors^. 
It is clarified that a nonmonotonic temperature dependence of the critical current appears in D/DN/D junctions 
due to the competition between the proximity effect and the formation of MARS. Concerning unconventional su- 
perconducting junctions with clean ferromagnet, Josephson effect is studied in Refr^i. However Josephson effect in 
unconventional superconducting junctions with diffusive ferromagnet has not yet been studied, although proximity 
^ ' effect in these junctions was studied recently^. Moreover detailed derivation of the boundary conditions is not given 
5^ ' and only D/DN/D junctions are considered in Ref.'^^. Since proximity effect and MARS strongly influence the den- 
sity of states'^^''^^''^'^, they should also crucially influence Josephson effect in diffusive unconventional superconducting 
junctions. 

The purpose of the present paper is to study the Josephson effect in various types of conventional and unconventional 
superconducting junctions with DN or DF interlayers. We flrst provide the technical details of derivation in Rcf.'^-. 
Then, solving the Usadel equations, we apply the general approach to study the influence of the exchange field, 
the proximity effect and the formation of MARS on Josephson current simultaneously. A number of peculiarities in 
the Josephson current are found depending on the pairing symmetry: a nonmonotonic temperature dependence in 
D/DN/D junctions, anomalous current-phase relation in P/DN/P junctions, second harmonics in the current-phase 
relation and their half-periodic oscillations as a function of the length of DF in D/DN/D and P/DN/P junctions. 



o 



2 



transitions to a tt state in D/DF/D, P/DF/P and S/DF/P junctions, a double peak structure in the temperature 
dependence of the critical current in D/DF/D junctions, and enhancement of Josephson current by the exchange field 
in S/DF/P junctions. 

II. FORMULATION 

We consider a junction consisting of unconventional superconductors (USCs) connected by a quasi-one-dimensional 
DN (or DF) with a resistance Rd and a length L much larger than the mean free path. The DN/USC interface 
located at a; = has the resistance while the DN/USC interface located at a; = L has the resistance Rb- We model 
infinitely narrow insulating barriers by the delta function U{x) = HS{x — L) + H'S(x). The resulting transparencies 
of the junctions T„ and are given by = 4cos^ 0/(4cos^ + Z^) and = 4cos^ 0/(4cos^ c/) + Z'^), where 
Z = 2H/vf and Z' = 2H' /vp are dimensionless constants, is the injection angle measured from the interface normal 
to the junction and vp is Fermi velocity. 

In order to study the Josephson effect in diffusive USC junctions, we first concentrate on the quasiclassical Keldysh- 
Nambu Green's function in DN defined by Gn{x). Its retarded part Rn{x) can be expressed as 

Rn{x) = cos ip sin Oti + sin ijj sin 9t2 + cos Ot^ , (1) 

with Pauli matrices in the electron-hole space, fi, f2, and f^. Since R^ix) obeys the Usadel equation, following 
equations are satisfied, 

D\—e-(^)^cos0sm9] + 2i(e+(-)h)sme^O, (2) 
ax'^ ox 

for majority (minority) spin with the diffusion constant D and exchange field h. The boundary condition of Gn{x) 
at DN/USC interface is given hy^l^ 



L - dGNjx) < > 
-[G^(x)^^]|..^_ = 



[Gn{x) 7J^]|x=L_ = ^_ , (4) 



/(0) = 2[Gi,B], (5) 
B = {-Ti[GuHZ^]+ HZ^H+ -TfGiHZ^H+Gi)-\Ti{l - HZ^) +TfGiHZ^H+) (6) 



with Gi = Gn{x = H± = (G2+±G2-)/2, and Ti = r/(2 - T + 2V]~T), where G2± is the asymptotic Green's 
function in USC as defined in our previous papersSI. The average over the various angles of injected particles at the 
interface is defined as 

/ir/2 (.7r/2 
d0 COS 0/(0)/ / d0T(0)cos0 (7) 
-■it/2 J-tt/2 

with /(0) / and r(0) = T. The resistance of the interface R\j is given by 

Rb — —TTT, ^ (8) 

n;2d0r(0)cos0 

with Sharvin resistance at the interface, Rq. Retarded components of Gi and G2± are given by i?i,2± where i?2± is 



expressed by R2± = g±T3 + /±T2 with g± = e/y _ /\2_ g^^^j _ ^ „ ^2^ where e denotes the quasiparticle 
energy measured from the Fermi energy. A_|_ (A_) is the effective pair potential felt by quasiparticles with an injection 
angle (tt — 0). We also denote i?p, Rm and In as retarded part of H- and I. 
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Now we discuss the boundary condition of the retarded part of Keldysh-Nambu Green's function at DN/USC 
interface. The left side of the boundary condition of Eq. (4) can be expressed as 



^-RNix) — RN{x)\j:^L = ;^[(^ — sin 6* cos 6* cos V' J fi+|^— cosV'- — sin 6* cos 6* sin J -r2 + — sin Ofs] 

(9) 

In the right side of Eq. (4), Ir can be expressed by using several spectral vectors: 

Ir = UTiidR ■ + T^)is2+ ~ S2- f[si X (S2+ + S2-)] ' T 

+ 2T1S1 • {S2+ X S2-)[SI X (S2+ X S2-)] ■ T 
+ 2TiSi ■ {S2+ - S2-)[SI X {S2+ ~ S2-)] ■ f 
+ Tl){\ - S2+ • S2-)[Si X (S2+ X S2-)] ■ T 
+ 2iTi(l - S2+ ■ S2-)[SI ■ {S2+ - S2-)SI - (S2+ - S2-)] ■ t}, (10) 
dR={l+ TI){S2+ X S2_) - 2riSi X {S2+ - S2-) - 2X^8^ ■ (^2+ X S2_)Si (11) 

with Ri = Si ■ T and i?2± — ■S2± • "t— • 

The spectral vectors Si and S2± are given by 

(sin 6 cos -0 \ / /± "^o*^ * 

sin6'sinV' , S2± = f±sm^ | (12) 
cos6' y \ 9± 

where denotes the phase of the USC. Then Ir is reduced to 



Ir ^2iT[{2^T) +T{fs sin 9 cos (V' - *) + gs cos 6*} - T/s sin 6* sin (V^ - *)] ^ 

X I [ —gs sin 9 sin ip + fs cos sin ^' — /s cos cos ^I'] f i 
+ [ — /s cos cos + gs sin cos ip ~ fs cos sin f2 

+ [/s sin 6* sin (-0 - + /s sin 9 cos (?A - *)] fa } , 

and hence we find the following form of the matrix current: 

Ii sin 9 sin -0 + -^2 cos 9 sin ^' — /a cos 9 cos 4* 
-^/?)=j( —^2 cos 6* cos ^ + /i sin cos — /a cos 6* sin 'J | • f, 
/2 sin 61 sin (-0 - vl/) + /g sin 6* cos (V' - 

J _ , 2T'm.9s \ J _ / '^T„ifs \ J _ / 2T„i,fs 



Finally the boundary conditions are given by 
Li?;, 5 



(13) 



A / ' ^ \ A / ' " \ ^ 
A ^{2- T,n) + T,n {fs sin 9 cos (^ - + .gs cos 0} - T„Js sin 9 sin (V^ - *), 



= + = + fs ^ '(f+9--f-g+) . . 

l + f+f-+9+g-''' 1 + /+/_+ 3+5- 1 + f+U + g+g-- ^ ^ 



= -/i sin 6* + /2 cos 6* cos (?A - *) - /a cos 6* sin {^p - *) , (15) 

^ sin 6*— = -/2 sin (V" - *) - /a cos (V' - ■ (16) 

Rd 8x 
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For the calculation of the thermodynamical quantities, we use Matsubara representation: e — > iuj. We parametrize 
the quasiclassical Green's functions G and F using function <&: 

G^ ^ ==cos6i, (17) 

= , ,^ \ ^'^G^^siuOe-'^, (18) 

Cj)* $* 

FI„ = = = ^^G^ = sin6'e*'^ (19) 



with Matsubara frequency uj. Then Usadel equation read; 



,38 



.2^Tc d / 2 d 



for majority (minority) spin with ^ = D/2tiTc and critical temperature Tc- The following relations are satisfied: 

sin cos = §^ (^-c. + , (21) 



2uj 

sine sin V = ^(^^-I'l^). (22) 



Then the boundary condition is expressed as 



Gui d ^ Rd I , 



UJ ax RhL \ (jj 



A /'^ \ A /'^ \ A 
A = 2 - T„ + Tra{gsG^ + Js [B cos * + G sin ^) - /s (G cos * - B sin ^)), 

S = |^(<i>. + $1J, C=^($.-$1J (23) 

at a:; = L. 

This boundary condition is quite general since with a proper choice of A±, it is applicable to any unconventional 
superconductor with 5*2 = in a time reversal symmetry conserving state. Here, Sz denotes the z-component of the 
total spin of a Cooper pair. For s-, d- and p-wave superconductors we choose A± = A(r). A(T) cos(2(/> =F 2a) and 
A(r) cos(0 =F ct) respectively. 

In the following we will calculate Josephson current using this boundary condition at a: = and x = L, where ip is 
the external phase difference across the junctions, and a and /3 denote the angles between the normal to the interface 
and the crystal axes of USCs for a; < and x > L respectively. It is important to note that the solution of the Usadel 
equation is invariant under the transformation a —a or /3 — s- — /3. This is clear by replacing (f) with — (/> in the 
angular averaging. 

Josephson current is given by the expression 

elR^^RTI^yGl/^ ^ \ 



with temperature T and R — R^ + Rb + R'l,. Below Ic denotes the critical current and we consider symmetric barriers 
with Rb = R'f^ and Z — Z' for simplicity. 
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III. RESULTS 



A. Junctions with DN 



Let us first focus on the junctfons with DN. Figure [T] shows the current-phase relation for T/Tc = 0.1, Rd/Rb = 0.1 
and i?T/i/A(0) = 0.2 in (a) s-wave, (b)(i-wave and (c)p-wave superconducting junctions with {a, 13) — (0,0). In 
s-wave junctions, the IR product is suppressed with the increase of Z because proximity effect is suppressed. In 
d-wave junctions, proximity effect and hence IR are enhanced with the increase of Z because of the cancellation of 
the positive and negative parts of pair potential in the angular averaging. As Z increases, the contribution from 
the positive part exceeds that from negative part and hence the cancellation becomes weak"^^. In p-wave junctions, 
the IR product is strongly enhanced with the increase of Z because of the formation of the resonant states. It is 
known that the proximity effect and MARS can coexist"^^. At {a, f3) — (0,0), proximity effect is mostly enhanced. 
As Z increases, the contribution of the MARS becomes remarkable and hence the proximity effect gets strongly 
enhanccd"^^. Consequently its magnitude is an order of magnitude larger than that in s-wave junctions. The results 
at lower temperatures are shown in Fig. [2] The IR product is enhanced with decreasing temperature because the 
proximity effect is enhanced. In this case the current-phase relation in p- wave junctions has the form close to sin(/3/2, 
in contrast to the standard sinusoidal relation. This is a peculiar property of the formation of the resonant states in 
p-wave junctions where constructive interference occurs near (p — 7r4^ 
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FIG. 1: (color online) Current-phase relation for T /Tc = 0.1, Rd/Rb = 0.1 and Exh/^i^) = 0.2. (a)s-wave junctions. 
(b)ci-wave junctions, (c)p-wave junctions. 

Next we study the Josephson effect for other misorientational angles. As a or /3 increase, the IR product is 
monotonically suppressed due to the suppression of the proximity effect as shown in Fig. [3l[a) and Fig. ^h). In 
d(p)-wave junctions, the first harmonics disappear and hence the IR product is proportional to — sin2</? at (a,/3) — 
(7r/4, 0) ((7r/2, 0)) as shown in Figs. [3] (c) and (d). This can be explained in the limiting case as follows. Near Tc, the 
IR product, which stems from the first harmonics, is proportional to cos 2a cos 2/3 in d-wave junctions because angular 
averaging gives < cos(2(/)— 2a) >oc cos2a.^ Thus the first harmonics disappear at a = 7r/4. Similar argument is also 
applicable to p-wave junctions. 

Let us discuss the results for the critical current. In Fig. 31 temperature dependence of the critical current is 
plotted for various Rd/Rb and i?Th/A(0) in (a) s-wave, (b) d-wave and (c) p-wave superconducting junctions with 
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FIG. 2: (color online) Current-phase relation for T /Tc = 0.02, Rd/Rt ~ 0.1 and Exh/^iO) = 0.2. (a)s-wave junctions. 
(b)d-wave junctions, (c)p-wave junctions. 

Z = 10 and (a,/3) — (0,0). As Rd/Rb and i?Tft/A(0) increase, Ic'R increases for all the junctions because proximity 
effect is enhanced. In p-wave junctions, the critical current is strongly enhanced at low temperatures compared to 
the s-wave and d-wave junctions. When the misorientational angles in d-wave junctions are changed, nonmonotonic 
temperature dependence appears as shown in Fig. [5] (a). This nonmonotonic behavior can be explained in terms of 
the competition between the proximity effect and the formation of MARS. It is known from the previous studies that 
for a = P = the proximity effect exists but MARS is absent at the interfaces. On the other hand, for a = /? = 7r/4, 
only MARS exists and the proximity effect is absen t*^^'^^ . In other cases, both the proximity effect and MARS are 
present. With the decrease of temperature, the formation of MARS strongly suppresses the proximity effect. This 
results in the suppression of the Josephson current at low temperatures. Therefore, a nonmonotonic temperature 
dependence appears when both the proximity effect and MARS coexist. 

The above statement can be confirmed by calculation of the dependence of anomalous Green's function F on 
Matsubara frequency ui as shown in Figs. (b) and [5] (c) at x — L/2 and ip — tt/2 for {a, f3) — (7r/8,0). At 
low temperature [T /Tc = 0.01) the magnitude of Imi^ is suppressed at low energy in contrast to the case of high 
temperature {T /Tc = 0.2 and 0.3). This result illustrates strong suppression of the proximity effect by the formation of 
MARS at low T, which leads to the nonmonotonic temperature dependence. Note that this nonmonotonic dependence 
can appear only for large Z when the role of MARS is essential 

It is interesting to study the junctions composed of superconductors with different symmetries. Here we study 
S/DN/D junctions with Z = 10, Rd/Rb = 1 and Eth/ ^i*^) = 0.1. We choose Tcd/Tcs = 5 in Fig. [SJa) and 
Tcd/Tcs = 10 in Fig. EUb) where Tcs{Tcd) denotes the critical temperature of the s-wave (d-wave) superconductors. 
In this case the nonmonotonic temperature dependence also occurs due to the competition as shown in Fig. [S] In S/I/D 
junctions, the nonmonotonic temperature dependence was observed experimentally in Ref.— . We can qualitatively 
explain these data by regarding the barrier as a diffusive normally conducting material. 

We study dependence of the critical current on barrier thickness L at various temperatures in (a) s-wave, (b)c?-wave 
and (c)p-wave superconducting junctions with Z — 10, Rd/Rb = 0.1 and (a, /?) — (0, 0) in Fig. [T] The IcR product is 
proportional to exp{—CL/S/) for large L/^ for all the junctions as shown in Fig. [71 Here C is a constant independent of 
L. As temperature is lowered, the magnitude of C is reduced. From our results, we also find the relation C cx T"^/^. 
The results for the junctions with other misorientational angles for T /Tc = 0.01, Z = 10 and Rd/Rb = 0.1 are shown 
in Fig. [H As a and /? increase, IcR is suppressed. However, C is independent of these values, which indicates that 
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FIG. 3: (color online) Current-phase relation for T/Tc = 0.01, Z = 10, RdjRh = 0.1 and £Th/A(0) = 0.2. (a) and (c)d-wave 
junctions, (b) and (d)p-wave junctions. 

MARS don't influence the effective coherence length f^jC. This is because the effective coherence length reflecting 
the penetration of Cooper pairs is determined by the Usadel equation and therefore is independent of MARS. 

B. Junctions with DF 

Here we consider junctions with DF. We will study three types of junctions: D/DF/D, P/DF/P and S/DF/P 
junctions. Figure [H shows current-phase relation in D/DF/D junctions for T jTc — 0.01, Z = 10, RdjRh = 1 and 
Erh/^iO) = 0.1. At (a,/3) = (0,0) where the MARS are absent, the exchange field causes a O-tt transition as 
predicted for s-wave junctions (see Fig. [3] (a)). Similarly, second harmonic changes its sign at {a, f3) = (7r/4, 0), where 
the proximity effect is absent at x = 0, as shown in Fig. [5] (b). Figure \TU\ displays temperature dependence of the 
critical current in D/DF/D junctions with Rd/Rb — 1 and Erh/^iO) = 0.1. At {a, /3) = (0,0), the exchange field 
causes a O-tt transition as shown in Fig. [TOKa). At (a, (3) = (tt/S, 0), the exchange field also causes a O-tt transition, and 
as a result, double peak structure appears for h/A{0) — 0.4 as shown in Fig. [TO] (b). The peak at lower temperature 
stems from the competition between proximity effect and MARS. The peak at higher temperature stems from the 
O-TT transition. With the decrease of Z, the magnitude of IcR is suppressed while the O-tt transition temperature is 
almost independent of Z (see Fig. [10] (c)). At (a,/3) = (tt/S, 0), the peak at lower temperature disappears for small 
Z as shown in Fig. [TO](d) because the existence of the insulating barrier is essential for the formation of MARS. 

The barrier thickness dependence of the critical current in D/DF/D junctions is plotted in Fig[TT]with T/Tc = 0.1, 
Z = 10, Rd/Rb = 1 and ETh/A{0) = 0.1. For h = 0, the IcR product has an exponential dependence on L. As 
h/A{0) increases, Ic oscillates as a function of L/£^. The period of the oscillation becomes shorter with increasing h as 
shown in Fig[TT](a). As a and /3 increase, IcR is suppressed while the period of the oscillations remains constant, that 
is, the period is independent of the MARS. The second harmonics have a shorter (almost half) oscillation period than 
that of the first harmonics, similar to the predictions for S/DF/S junctions^^'"^^ (see the result for (a,/3) = (7r/4,0) in 
Fig[II](b)). 

Next we consider the P/DF/P junctions. Current-phase relation in P/DF/P junctions for T/Tc = 0.01, Z — 10, 
Rd/Rb = 1 and Eth/ AiQ) = 0.1 is plotted in Fig[TO] With increasing h, the dependence of IR changes from sin(p/2 
to sin 293 and finally to — sintp at {a, /3) = (0,0) as shown in Fig[T2](a). The phase dependences originate from the 
formation of the resonant states, the disappearance of the first harmonics at the O-tt transition and the emergence of 
the TT-junctions, respectively. At (a, f3) = {tt/2, 0) where the MARS are absent at x = 0, the second harmonics change 
the sign with the increase of h as shown in Fig[TO] (b). The critical current as a function of T is shown in Fig[T3] (a). 
The O-TT transition occurs due to the exchange field. Similarly, as h/A{0) increases, IcR oscillates as a function of 
L/S^. The period of the oscillation becomes short with increasing h as shown in FigfTOI (b). 

Finally we study S/DF/P junctions for Z = 10, Rd/Rb = 1, Erh/ A(Q) = 0.1 and /3 = 0. Current-phase relation at 
T/Tc = 0.01 has the form of — sin 2{p for h = Q due to the difference of the parities of two superconductors^ as shown 
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FIG. 4: (color online) Temperature dependence of the critical current, (a)s-wave junctions, (b)d-wave junctions, (c)p-wave 
junctions. a..Ra/Rb = 2 and Erh/MO) = 1- h.Ra/Rt = 0.5 and STh/A(0) = 1. c.Ra/Rb = 2 and Eth/MO) = 0.1. 
A.Rd/Rb = 0.5 and Erh/HQ) = 0.1. 

in Fig |14l (a). As h increases, the shape of IR transforms from — sin2(/9 to cosLp since the first harmonics recover by 
breaking the symmetry between up- and down- spins. Temperature dependence of the critical current is plotted in 
Fig ll4l (b). The magnitude of IcR is enhanced by the increase of h due to the recovery of the first harmonics, in 
contrast to junctions between superconductors with equal parities. 

IV. CONCLUSIONS 

In this paper, we studied the Josephson effect in junctions between unconventional superconductors with diffusive 
barriers. The Usadel equations in the barrier region were solved with the generalized boundary conditions applicable 
to the unconventional superconductors at the interfaces. Applying these boundary conditions, we calculated the 
Josephson current in various types of junctions: S/DN/S, D/DN/D, P/DN/P, S/DN/D, D/DF/D, P/DF/P and 
S/DF/P junctions. Our main conclusions can be summarized as follows. 

1. The dependences of Josephson current on the interface barrier strength Z are different for S/DN/S, D/DN/D, 
and P/DN/P junctions. Josephson current is suppressed by the increase of Z in S/DN/S junctions while it is enhanced 
by the increase of Z in D/DN/D and P/DN/P junctions. In D/DN/D and P/DN/P junctions, proximity effect is 
enhanced by the increase of Z due to the cancellation of the positive and negative parts of pair potential in the 
angular averaging and the coexistence of MARS and proximity effect respectively. The coexistence also induces 
anomalous current-phase relation in P/DN/P junctions. When proximity effect is absent at one interface, the second 
harmonics dominate in D/DN/D and P/DN/P junctions. The competition between MARS and proximity effect 
causes a nonmonotonic temperature dependence of the critical current in D/DN/D junctions. Similar dependence can 
be seen in S/DN/D junctions. 

2. In S/DN/S, D/DN/D and P/DN/P junctions, the critical current has an exponential dependence on the length 
of the DN. The prefactor of the length is independent of the MARS. 
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FIG. 5: (color online) (a) Temperature dependence of the critical current for (a, /3) — (tt/S, 0). a..Rd/Rb ~ 2 and Eth/ ^ifi) ~ 1- 
h.Ra/Rb = 0.5 and £;Th/A(0) = 1. c.Ra/Rb = 2 and ETh/A{0) = 0.1. A.Rd/Rb = 0.5 and £Th/A(0) = 0.1. (b) real and (c) 
imaginary parts of anomalous Green's functions F with Rd/Rb = 2 and Erh/^iO) = 1. 



3. In D/DF/D, P/DF/P and S/DF/P junctions, the tt -state can be realized. A double peak structure in tem- 
perature dependence of the critical current occurs in D/DF/D junctions due to O-tt transition and the competition 
between MARS and proximity effect. In S/DF/P junctions, the Josephson current can be enhanced by the exchange 
field, in contrast to other types of junctions, due to the recovery of the first harmonics. 

4. In D/DF/D and P/DF/P junctions, the critical current has an oscillatory behavior as a function of the length 
of the DF. The period of the oscillation becomes short for large exchange field while it is independent of the MARS. 
The second harmonics show almost half periodicity compared to the first harmonics. 
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FIG. 8: (color online) Length dependence of the critical current with T jTc ~ 0.01, Z = 10 and Rd/Rb = 0.1. (a)d-wave 
junctions, (b)p-wave junctions. 




FIG. 9: (color online) Current-phase relation in D/DF/D junctions for T/Tc = 0.01, Z = 10, Rd/Rb = 1 and ETh/A{0) = 0.1. 
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FIG. 10: (color online) Temperature dependence of the critical current in D/DF/D junctions. Z = 10 in (a) and (b). /i/A(0) = 
0.4 in (c) and (d). 




FIG. 11: (color online) Length dependence of the critical current in D/DF/D junctions with T /Tc = 0.1, Z = 10, Rd/Rb = 1 
and Erh/ A{0) = 0.1. We choose ^^(O) = 4 in (b). 
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FIG. 12: (color online) Current-phase relation in P/DF/P junctions for T/Tc = 0.01, Z = 10, RdlRb = 1 and _ETh/A(0) = 0.1. 




FIG. 13: (color online) (a) temperature and (b) length dependence of the critical current in P/DF/P junctions. Z 
Rd/Rb = 1, and (q,/3) = (0,0). We choose -BTh/A(0) = 0.1 in (a) and T/Tc = 0.1 in (b). 



= 10, 
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FIG. 14: (color online) (a) Current-phase relation at T /Tc ~ 0.01 and (b) temperature dependence of the critical current in 
S/DF/P junctions for Z = 10, Rd/Rb = 1, £Th/A(0) = 0.1 and /3 = 0. 



